home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / chetri.z / chetri
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCHHHHEEEETTTTRRRRIIII((((3333FFFF))))                                                          CCCCHHHHEEEETTTTRRRRIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CHETRI - compute the inverse of a complex Hermitian indefinite matrix A
  10.      using the factorization A = U*D*U**H or A = L*D*L**H computed by CHETRF
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CHETRI( UPLO, N, A, LDA, IPIV, WORK, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, LDA, N
  18.  
  19.          INTEGER        IPIV( * )
  20.  
  21.          COMPLEX        A( LDA, * ), WORK( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      CHETRI computes the inverse of a complex Hermitian indefinite matrix A
  25.      using the factorization A = U*D*U**H or A = L*D*L**H computed by CHETRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              Specifies whether the details of the factorization are stored as
  31.              an upper or lower triangular matrix.  = 'U':  Upper triangular,
  32.              form is A = U*D*U**H;
  33.              = 'L':  Lower triangular, form is A = L*D*L**H.
  34.  
  35.      N       (input) INTEGER
  36.              The order of the matrix A.  N >= 0.
  37.  
  38.      A       (input/output) COMPLEX array, dimension (LDA,N)
  39.              On entry, the block diagonal matrix D and the multipliers used to
  40.              obtain the factor U or L as computed by CHETRF.
  41.  
  42.              On exit, if INFO = 0, the (Hermitian) inverse of the original
  43.              matrix.  If UPLO = 'U', the upper triangular part of the inverse
  44.              is formed and the part of A below the diagonal is not referenced;
  45.              if UPLO = 'L' the lower triangular part of the inverse is formed
  46.              and the part of A above the diagonal is not referenced.
  47.  
  48.      LDA     (input) INTEGER
  49.              The leading dimension of the array A.  LDA >= max(1,N).
  50.  
  51.      IPIV    (input) INTEGER array, dimension (N)
  52.              Details of the interchanges and the block structure of D as
  53.              determined by CHETRF.
  54.  
  55.      WORK    (workspace) COMPLEX array, dimension (N)
  56.  
  57.      INFO    (output) INTEGER
  58.              = 0: successful exit
  59.              < 0: if INFO = -i, the i-th argument had an illegal value
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCHHHHEEEETTTTRRRRIIII((((3333FFFF))))                                                          CCCCHHHHEEEETTTTRRRRIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              > 0: if INFO = i, D(i,i) = 0; the matrix is singular and its
  75.              inverse could not be computed.
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.